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f^ ■ Abstract. The Uhrig dynamical decoupling sequence achieves high-order decoupling 

^^ , of a single system qubit from its dephasing bath through the use of bang-bang 

Pauli pulses at appropriately timed intervals. High-order decoupling of single and 
i-Q . multiple qubit systems from baths causing both dephasing and relaxation can also 

be achieved through the nested application of Uhrig sequences, again using single- 
^ ' qubit Pauli pulses. For the 3-qubit decoherence free subsystem (DFS) and related 

C^ . subsystem encodings, Pauli pulses are not naturally available operations; instead, 

^ ' exchange interactions provide all required encoded operations. Here we demonstrate 

that exchange interactions alone can achieve high-order decoupling against general 
noise in the 3-qubit DFS. We present decoupling sequences for a 3-qubit DFS coupled 
to classical and quantum baths and evaluate the performance of the sequences through 
numerical simulations. 



PACS numbers: 03.67.Pp, 03.67.Lx, 03.65.Yz, 76.60.Lz 

1. Introduction 

Dynamical decoupling (DD) pulse sequences have had a long history beginning with the 
Hahn spin echo [1] and the Carr-Purcell-Meiboom-Gill sequences [2], |3] and continuing 
to the present. Recently, Uhrig has developed a DD sequence that, by varying the time 
intervals between pulses, is able to decouple transverse dephasing to order n in the 
system-bath coupling strength with n + 1 pulse intervals [H [5] . The Uhrig decoupling 
sequence is in fact universal, decouphng both classical and quantum baths that cause 
either transverse dephasing or longitudinal relaxation [6]. Several generalizations of the 
Uhrig sequence, and its use of non-uniform pulse intervals, have now been made. These 
generalizations allow high-order decoupling of a single qubit from dephasing noise baths 
with different noise power spectra [318], decoupling of a single qubit from baths causing 
simultaneous dephasing and relaxation [9] , and decoupling of multi-qubit systems from 
general baths [TO]. The Uhrig decoupling scheme and its generalizations all require 
single qubit Pauli pulses for implementation. 

In semiconductor quantum dot systems, however, single qubit operations are not 
easily implemented, in contrast with the two qubit exchange interaction. While the 
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exchange interaction can be performed in sub-nanosecond time scales [TT], single qubit 
rotations may be two orders of magnitude slower or more p!2| [T3] , and can be technically 
more demanding. The difference in requirements between single and two qubit gates has 
led to the development of encodings that use the exchange interaction alone [HI [151 in] • 
The smallest such encoding is the 3-qubit decoherence free subsystem (DFS), for which 
explicit exchange gate sequences for encoded universal computation are given in pTlflS] . 
Creation of a DFS itself may also be performed using exchange pulses alone [19] . While 
encoded computation can be performed with exchange gates alone, the use of the DD 
pulse sequences described above would require single qubit Pauli operations. Here we 
demonstrate that exchange gates alone suffice for high-order decoupling of the DFS- 
encoded information from general baths. 

Our new exchange-only DD sequences explicitly take advantage of the decoherence 
free properties of DFS encodings. Because the 3-qubit DFS is decoherence free with 
respect to any global interaction, decoupling from a decohering bath can be achieved 
by globalizing any local interactions to high order. In contrast to standard Pauli- 
based decoupling schemes, our aim is not to cancel system-bath coupling terms, but 
to globalize them, so that the effective Hamiltonian created by the DD scheme causes 
only the gauge qubit and bath to evolve. This idea applies to any subsystem encoding: 
decoupling need only preserve the subsystem of interest while gauge subsystems can 
evolve arbitrarily. An alternative method for low-order leakage elimination using the 
simultaneous operation of multiple exchange gates has been given in [20l[2T| . 

The paper is organized as follows. In section |2] we briefly review aspects of the Uhrig 
dynamical decoupling sequence and the 3-qubit decoherence free subsystem. Section 
E] describes how exchange-only decoupling is achieved for a 3-qubit DFS subject to 
dephasing from a classical noise bath. Section H] gives the analogous exposition for 
decoupling from a quantum bath. Simulations of the decoupling sequences for a DFS 
qubit coupled to classical and quantum baths are presented in section [5l We conclude 
in section [61 In the following we use the acronym DFS to refer to decoherence free 
subsystem (rather than subspace). We assume "bang-bang" exchange pulses perform 
the decoupling, i.e., the pulses perform a perfect finite operation using infinite power in 
infinitesimal time. 

2. Background 

2.1. Uhrig dynamical decoupling 

The Uhrig dynamical decoupling (UDD) sequence decouples a single qubit from a 
classical dephasing bath [H [5] . The Hamiltonian for such a system is given by 

H = ZB{t), (1) 

where Z is the Pauli Z operator and B{t) is a time-dependent real- valued function. 
Subjected to a sequence of vr pulses about the x-axis, the Hamiltonian in the toggling 
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frame becomes 

H = fit)ZB{t), (2) 

where f{t) is the UDD "switching function" and takes values of ±1. The switching 
times, when / switches between 1 and -1 (or vice versa), give the instances when vr 
pulses are apphed. The propagator corresponding to free evolution interspersed with vr 
pulses is given by 

t/(r)=exp(-iZ^(T)), (3) 

where T is the total time of the DD sequence, and the accumulated phase 6 from the 
classical noise bath over time T is 

e{T) = [ dsf{s)B{s). (4) 

Jo 

For order n decoupling, we require that 

= / dsfis)sP (5) 

for all < p < n, which leaves only terms of order T"""*"^ and higher in the accumulated 
phase. In other words, the DD sequence has resulted in an effective Hamiltonian with 
system-bath coupling only at order T""*"^ and higher. Uhrig showed that the switching 
times for order n decoupling are given by the simple formula 

t,=Tsm^( /"" , V (6) 

' \2{n + l)J ' ^ ^ 

for j = 1, . . . ,n. It was subsequently shown that the UDD sequence is in fact universal 
[6], i.e., UDD decouples dephasing due to both classical and quantum baths. 

The propagator ([3]) evolves an initial state \ipi) to the final state \ijjf) = U(T) \ipi). 
The fidelity of memory preservation due to the decoupled evolution is then 

F{T) = |(^/|^.)P = cos'9{T) + {ij,\Z\i:,)'sm'e{T). (7) 

For an initial state on the y-axis of the Bloch sphere (ipilZlipi) = 0, and the performance 
of the DD sequence can be described in terms of the decoherence function W{T) = 
{Y{T)) [22], which is the expectation value of the Pauli Y operator. When {ifji\Z\tlJi) = 0, 
the fidelity and the decoherence function are related through F{T) = 1(1 + W{T)). For 
a noise bath with Gaussian statistics, the ensemble averaged decoherence function is 

W{T) = e-^(^\ (8) 

X(T) = r^5(a;)|/»(u;T)|l (9) 

Jo ^TT 

Here S{u) is the bath's noise spectral density and /'•"^(wT) is the Fourier transform of 
the UDD switching function f^"'\t) of order n: 

2 
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Since {tj} are the UDD switching times, {tk-i+tk)/2 is the midpoint of the k^^ interval, 
where /'•"^ is (—1)'^"^. The expression in fITOl) is the "filter function"; together with ([8]) 
and ([9]) the filter function gives an interpretation of the UDD sequence as a solution 
to a filter design problem [23]. For the UDD sequence, the filter function gives a high- 
pass filter, suitable for decoupling noise spectra with a sharp high-frequency cutoff. 
The switching times can be adjusted to decouple noise spectra with other frequency 
characteristics [71 [8] . 



2.2. The 3-qubit decoherence free subsystem 

The general theory of decoherence free subspaces and subsystems and encoded 
universality is given in [15]. Initialization, measurement, and universal computation 
in the 3-qubit DFS are described in [TTJIIH]. In this subsection we give a brief summary 
of the 3-qubit DFS to fix notation and to provide the physical motivation behind the 
pulse sequence designs in the following sections. 

The states of three spin-^ qubits can be described by the quantum numbers of 
three commuting operators S, 5*1,2, and Sz- S is the total spin of all three qubits and 
distinguishes valid and leaked states in the DFS encoding. 5*1,2 is the total spin of the 
first two of the three qubits and determines the encoded qubit state. Sz is the total 
2;-axis spin of the three qubits and gives the gauge state. Exphcit definitions of these 
operators are given in [18]. The states of the eight dimensional Hilbert space of the 
three qubits are explicitly: 
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where states on the left in (TTT]) - (TTS1) are assigned labels 1-8, states in the middle are 
described in the angular momentum basis |5, 5i,2, Sz), and states on the right are written 
in terms of the standard computational basis. 
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States |1)-|4) in (TTT]) - (IT^ span the valid subspace of the 3-qubit DFS. States |1) 
and |2) are encoded states, with gauge states +\ and — |, respectively; states |3) 
and |4) are encoded 1 states, with gauge states -\-\ and — |, respectively. In the valid 
subspace, the encoded quantum number 5*1^2 and the gauge quantum number S^ may be 
considered as the computational basis states of two two-state subsystems — two effective 
qubits. The states have been ordered so that the first subsystem effective qubit gives the 
encoded state in the DFS and the second subsystem effective qubit gives the gauge state. 
Valid DFS states must not only be in the valid subspace, but additionally be factorizable 
between the encoded and gauge quantum numbers. We also define a projector onto the 
valid subspace 

n^|l)(l| + |2)(2| + |3)(3| + |4)(4|, (19) 

for use in section [31 

Encoded operations are performed with exchange gates between pairs of constituent 
physical qubits. Because exchange commutes with the total angular momentum 
operators S and Sz-, exchange can only change the encoded DFS quantum number while 
leaving the gauge and leaked quantum numbers unaffected. Similarly, any interaction 
comprised of total spin component operators can only change the gauge quantum number 
Sz and leaves the encoded S'1,2 and leaked S quantum numbers unchanged. The encoded 
qubit is thus decoherence free with respect to any global interaction. For a given system- 
bath coupling of the DFS constituent physical qubits, our goal is to design a pulse 
sequence resulting in an effective Hamiltonian whose system-bath interactions to order 
n contain only total spin component Sx, Sy, or Sz operators. 

3. Decoupling the 3-qubit DFS from classical phase noise 

Consider the following Hamiltonian coupling a 3-qubit DFS to a classical constant 
dephasing noise bath: 

H = Hi = ZiBi + ^2^2 + ^3^3, (20) 

where Zj gives the Pauli Z operator for physical qubit j, and Bj is a constant (in time) 
real number giving the bath strength at qubit j. Free evolution corresponding to (!20l) 
is Uf{t) = exp(— lift). We define the pulse P as 

P = P2,3-Pl,2 (21) 

where Pij is a full swap operation between qubits i and j. The Hamiltonian in (l20l) 
conjugated singly and doubly by the pulse P gives the permuted Hamiltonians 

H2 = P-\H^.P = Z2B^ + ZsB2 + Z,B■,, (22) 

H:i = P-\H2.P = Z3B1 + Z1B2 + Z2BS. (23) 

Conjugation by P causes the physical qubits to dephase under the influence of the 
different local noise baths. By spending equal time in each of the three permutations 
in (l20l) . (!22l) . and (!23|) — the even permutations (alternating group ^3) of the symmetric 
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group 5*3 on three elements — we generate an effective global system-bath interaction. 
The propagator corresponding to equal time evolution in each permutation is given by 

U = expi-iHsT). exp{-iH2T).exp{-iHiT) (24) 

= exp(-i(Zi + Z2 + Zs)iB, + B2 + B,)t) (25) 

= P.Uf{T).P.Uf{T).P.Uf{T). (26) 

The effective Hamiltonian in ( !25l) consists of the total spin-z operator alone and evolves 
only the gauge qubit. We have thus decoupled the encoded subsystem state from the 
dephasing bath to first order. The pulse sequence for this first order exchange-only 
scheme is shown explicitly in ( l26l) and involves free evolution and full swap (exchange) 
operations alone. Of course the bath functions are not generally constant, and in the 
following we describe how this "globalization" of the Hamiltonian can be accomplished 
to high order for general time-varying bath functions. We will see that free evolution 
for non-uniform time intervals, interspersed with full swap operations, is sufficient for 
high order decoupling. 

For local classical dephasing baths with arbitrary time dependence the Hamiltonian 
is 

H = ZiBiit) + Z2B2{t) + Z,Bs{t). (27) 

Since the effect of conjugation by P or P^^ is to permute the local bath seen by the 
constituent qubits, the Hamiltonian in the toggling frame (cf. ([2])) is 

H = ZiBa^^t){t) + Z2B^,(^t){t) + Z3Ba,^t){t), (28) 

where ctjit) = 1,2, or 3 and identifies which local bath constituent qubit j is 
experiencing. Restriction to ^3 permutations completely specifies 02 (^) and 03 (t) in 
terms of ai{t): 02 (t) = «i(t) + l(mod 3) and 0:3 (t) = ai(t) + 2 (mod 3), where the 
modulus is taken with offset 1. For notational convenience we use all three aj(t)'s in the 
following. The three types of Hamiltonians Hi (120|) , H2 (!22l) , and H3 (123!) (generalized 
to time-dependent local baths) have ai values of 1, 3, and 2, respectively. 

The Q;j(t)'s change values when P or P~^ pulses are applied. Between pulses the 
cnjit) values are constant. The instances at which the pulses are applied are the switching 
times for exchange-only decoupling. Unlike the UDD situation, both the switching times 
and the Hamiltonian sequence (or ai{t)) must be determined. The UDD case involves 
only two Hamiltonian types and one pulse type: the Pauli X pulse toggles the evolution 
back and forth between the two Hamiltonian types. In the DFS case, we may choose to 
apply either a P or P~^ pulse, i.e., after evolution under one of the Hamiltonians (l20l) . 
(!22|) , or (!23|) , the next evolution interval can be chosen from either of the remaining two 
Hamiltonians. 

The propagator associated with the modulated Hamiltonian fl28l) is 

U{T) = exp i-iZMT)) exp (-iZ2^2(T)) exp {-iZ-MT)) , (29) 

where Oj{T) gives the phase accumulated by constituent qubit j in time T, 



e,iT) = / dsB^^^s)is). (30) 
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If all the 9j{T) values are equal, the propagator will contain only a global interaction 
term generated by the total spin 5*^, which drives gauge evolution only. Computing the 
fidelity of memory preservation shows explicitly that setting all the Oj{T) equal leaves 
the encoded subsystem unchanged. Given an initial DFS state \ipe) j-^g), where |'?/'e) 
gives the encoded subsystem state and \ipg) gives the gauge subsystem state, the fidelity 
of memory preservation under the evolution in fl29|) is 

i^m = EKH(^e|^f/(T)^|v^e)|V^,)|^ (31) 

The sum over /x is a partial trace over the gauge subsystem, since we are concerned only 
with encoded subsystem preservation. 11 is the projector into the valid subspace defined 
in ( IT9|) . Substituting the propagator (129|) into the fidelity ( 13T|) we find 

F = Co + ci cos2(^2(T) - e;{T)) + C2Cos2{e^{T) - e,{T)) 

+ C3cos2(^i(r)-02(r)), (32) 

where the c/s depend on the initial encoded state \ipe) (and not the initial gauge state 
IV^g)) and are given in the Appendix; the c/s sum to 1. 

Equation f p2|) shows that setting all the accumulated phases equal preserves the 
encoded subsystem fidelity. From (l30l) the accumulated phase difference 6*1 (T) — 6*2 (T) 
is 

e,{T)-e2{T)= [ ds(5„,(,)(s)-5,,(,)(s)) (33) 

Jo 

= [ ds{h{s)B,{s) + Ms)B,{s) + f,{s)B,{s)). (34) 

Jo 

The switching functions fj{t) depend on the Hamiltonian type in each evolution interval. 

For evolution under Hi, qubit 1 sees bath Bi and qubit 2 sees bath B2 so ai = 1 and 

0^2 = 2. The phase difference between qubit 1 and qubit 2 gives a positive sign on i?i, a 

negative sign on B2, and no contribution from B3 so that /i = 1, /2 = — 1 and /a = 0. 

The switching functions and aj values for the other Hamiltonian types may be similarly 

determined and are given in table [H The two remaining phase differences are 

e2(T)-^3(T)= / dsiMs)Biis) + fiis)B2is) + Ms)B,is)), (35) 

Jo 

e,{T)-ei{T)= [ dsiMs)Biis) + his)B2is) + Ms)B,is)). (36) 

Jo 

The same switching functions fj{t) appear here as in (IMj) . 

Decoupling of the DFS from arbitrary time-dependent classical baths can be 

achieved by setting each term in ( 1M1) . ( l35l) . and ( l36l) to zero, order-by-order. For order 

n decoupling this results in the constraints 

0= / dsMs)sP, (37) 

Jo 

0= [ dshis)sP, (38) 
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Table 1. Bath identification functions aj and switcliing functions fk- 



Hamiltonian type 


oil 


a2 


as 


/i 


/2 


/3 


Hi 


1 


2 


3 


1 


-1 





H2 


3 


1 


2 


-1 





1 


H3 


2 


3 


1 





1 


-1 



ioT < p < n and a suitably chosen sequence of Hamiltonian types. Because 
/s = — (/i + /2), satisfying f p7|) and fl55]) automatically constrains /^ ds/3(s)s^ to be 
zero. Equations (137|) and ( 138|) double the number of UDD constraints ([5]), and use two 
different switching functions that take values 1,-1, and 0. 

We have found numerical solutions to the constraint equations ( 137|) and fl38|l to 64**^ 
order using a basic period four sequence of Hamiltonians {Hi,H2,H^,H2\. Switching 
from Hi to H2 and from H2 to iJs is accomplished with a P pulse; P~^ brings H^ 
to H2 and H2 to Hi. The basic period four sequences for the switching functions are 
/i = {1, —1, 0, —1} and /2 = { — 1, 0, 1, 0}. For order n decoupling the basic period four 
sequence is repeated \{2n + l)/4] times, and the first 2n + 1 Hamiltonians are taken. 
The sequence of pulses P and P~^ is determined from the progression of Hamiltonian 
types. Odd decoupling orders require a final P pulse at the end of the sequence 
(in general, the product of the pulses over the whole evolution interval must be the 
identity). For example, third order decoupling repeats the sequence twice, and takes 
the first seven entries {Hi,H2,H3,H2, Hi,H2,H^}. The corresponding pulse sequence 
is {J-(ri),P, J-(r2),P,^(r3),P-S^(r4),P-S^(r5),P,^(r6),P,^(r7),P}, where ^(r) 
is a free evolution interval of length r. The time intervals Tj are given by the difference 
tj — tj-i between successive switching times. Figure [U shows the ^3 DFS DD switching 
times (filled circles • ) up to 10*^ order satisfying f l37|) and fl38|) . The UDD switching 
times (open circles O ) are also shown; a pair of ^3 DFS DD switching times brackets 
each of the UDD switching times. Though this bracketing structure is suggestive, we 
have been unable to determine a simple analytical formula for the DFS DD switching 
times. Numerical values for these switching times are given in table El 

Filter functions corresponding to the Fourier transform of the two order n A^ 
switching functions fi (t) and f2v) can be computed analogously to (fTOl) . Figure 
|2] shows the filter functions corresponding to the UDD switching function and the A^ 
DFS DD switching functions for various orders of decoupling. Since all the switching 
functions of a given order are designed to integrate the same number of monomials to 
zero, they all exhibit the same low frequency filtering behavior. 

We have thus far restricted ourselves to decoupling over even permutations ^3 only. 
The previous considerations can be generalized to decoupling over all permutations of 
5*3. Defining the odd permutation Hamiltonian types as 

H^ = Z2B1 + Z1B2 + Z3P3, (39) 

H^ = Z3B1 + Z2B2 + Z1B3, (40) 
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Figure 1. Switching times for orders n = 1 to n = 10. UDD switclring times are open 
circles O . Even permutation A3 DFS DD switcliing times are filled circles • only. 
Switching times for ^3 DFS DD over all six 6*3 permutations include open circles, filled 
circles, and stars *. The switching times are symmetric about i = 5. 
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Figure 2. Filter functions for various decoupling orders. Dashed lines are the UDD 
filter functions. Black and gray lines are the A^ DFS DD filter functions corresponding 
to /{" (t) and /2 (t) , respectively. The sets of lines are labeled by their decoupling 
order n. n = is free evolution. 
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Table 2. A-^ DFS DD switching times corresponding to the fihed circles in figure [T] 

Switching times hsted are for the interval [0, 1]. Only switching times less than i are 
given explicitly. Switching times greater than ^ are obtained by reflecting the given 
values about i = i. 



order n A^ switching times order n A3 switching times 



1 



0.3333333333333333 



0.1666666666666667 
0.3333333333333333 

0.0930802599812912 
0.2041913710924023 
0A4UUUUUUU 

0.0611678063574247 
0.1320291453900112 
0.2986958120566778 
0.3945011396907580 

0.0422244245173296 
0.0940587956886883 
0.2172228408817372 
0.2838895075484039 
0.4518343711713587 

0.0313685011617312 
0.0691609286752199 
0.1617103538537611 
0.2161866929592387 
0.3514848584641742 
0.4258827585118745 

0.0239219438795333 
0.0535688803938237 
0.1262566342290569 
0.1675244212375237 
0.2761133079137736 
0.3417044666375784 
0.4698392155798953 



8 0.0190156712850090 
0.0422945303794296 
0.1002297726086257 
0.1346268067223472 
0.2239571558152790 
0.2763447583052162 
0.3850761827426867 
0.4416793537112741 

9 0.0153608717513108 
0.0344809416081787 
0.0820319124861268 
0.1096019013513601 
0.1835330371665574 
0.2291713148467980 
0.3223652733904291 
0.3688693585992699 
0.4721657549445159 

10 0.0127428989292003 
0.0284688256262034 
0.0679240161384205 
0.0914464121824144 
0.1538757061482468 
0.1916101903303824 
0.2714006848897883 
0.3135792550438800 
0.4050737288155140 
0.4525790184049564 



H. 



Z,B, + ^3^2 + ^2^3, (41) 

and computing the accumulated phase difference integrals (cf. flM|) - fl35|) ) 
give rise to five sets of constraint equations (cf. ( 1371) and ( l38l) ) and five 
independent S3 switching functions. We find that a basic period 10 sequence of 
{Hi, H4, H2, H^, H3, Hq, H3, H^, H2, H4} allows order n decoupling with 5n + 1 evolution 
intervals. The period 10 sequence alternates between even and odd permutations, with 
the transformation between successive Hamiltonians accomplished with a single swap 
gate on qubits 1 and 2, or qubits 2 and 3. The corresponding basic sequence of swap 
gates is {-Pi,2, A.a, -^1,2, -^2,3, -^1,2, -Pi,2, -^2,3, -^1,2, -^2,3}- As before, the final pulse in the 
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sequence is chosen so that the pulses multiply to the identity. We can again solve 
the constraint equations numerically, and the switching times for 5*3 DFS DD are also 
displayed in figure [H (open circles O , filled circles • , and stars -k). The switching times 
for 5*3 DFS DD are not unique. Because decoupling can be accomplished with even 
permutations alone, or with odd permutations alone, the relative weighting between 
even permutation decoupling and odd permutation decoupling is unconstrained. The 
switching times for 5*3 DFS DD shown in figure [1] correspond to equal total time in even 
and odd permutations. Filter functions corresponding to the five S3 switching functions 
may be computed according to (fTOj) . They have the same low- frequency roll-off behavior 
as the UDD filter function. 

4. Decoupling the 3-qubit DFS from a quantum bath 

A DFS qubit coupled to local quantum baths is described by the Hamiltonian 

H = Bo + Si.Bi + S2.B2 + S3.B3. (42) 

Here Bq is a pure bath operator describing interactions within the bath alone, and 
Bj = {Bj^x, Bj^y, Bj^z) is a vector of bath operators coupled to constituent qubit j 
through its Pauh operators Sj = {Xj,Yj, Zj). As in section |3] full swaps are used to 
permute the qubits through the different bath operators. We define the quantum bath 
Hamiltonians analogously to (EO]), ^-^, and (I39l)-(HT1): 

H^ = Bo + SrA + 5^2.^2 + ^3.^3, (43) 

H2 = Bo + S2.Bi + Ss.B2 + S^.B3, (44) 

H3 = Bo + SsA + S,.B2 + S2.B3, (45) 

H, = B, + S2A + 5^1.^2 + ^3.^3, (46) 

H, = Bo + SsA + S2.B2 + Si.Bs, (47) 

He = Bo + Si.Bi + Ss.B2 + ^2.^3. (48) 

The first and second order A3 and 5*3 DFS DD sequences for decoupling classical 
phase noise also decouple the DFS qubit from quantum baths. Consider the first order 
A3 evolution sequence 

P. exp{-iHT).P. exp(-i/fr).P. exp(-iifr) (49) 

= exp(— ii73r). exp(— iiJ2T)- exp(— ii^ir) (50) 

^ 1 - ir (sBo + StotAot) , (51) 

where 

'S'tot = 5*1 + 5*2 + 5*3, (52) 

5tot = B^ + B2 + B3. (53) 

The first order effective Hamiltonian contains only coupling to the total spin component 
operators of all three qubits. Similarly, the second order ^3 evolution sequence is 

exp{-iHr).P'\ exp(-ii7r).P"\ exp(-2iifr).P. exp(-iifr).P. exp(-iifr) (54) 
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Table 3. 6*3 DFS DD switching times corresponding to the stars in figure [TJ Again, 
only switching times less than ^ are given explicitly (see the caption to table [2]) . Full 
53 switching times include the values below, the values in table [21 and the UDD values 

m_ 

order n S3 switching times order n S3 switching times 



1 



0.1666666666666667 



0.0833333333333333 
0.4166666666666667 



0.0441757320558095 
0.2663979542780318 
0.: 



0.0292438385042891 
0.1706622892054447 
0.2539956225387781 
0.4459105051709558 



0.0198486448526978 
0.1234090460471676 
0.1855655208780249 
0.3188988542113582 
0.4035604011944698 

0.0148169093658703 
0.0902375649702305 
0.1365285435153074 
0.2455460335064556 
0.3140624574739186 
0.4629576452117436 

0.0112075501170748 
0.0704161635276064 
0.1069644388480345 
0.1894875256294779 
0.2440254284265229 
0.3752930054921819 
0.4396659439243005 



8 0.0089377851765520 
0.0553969487226799 
0.0843798400197465 
0.1531864662990289 
0.1981928332289538 
0.3030094142447375 
0.3573598562810042 
0.4706108187731435 

9 0.0071855182206674 
0.0453635718581324 
0.0692317973728011 
0.1243583767250380 
0.1614266670585074 
0.2527437110834672 
0.2994843851270477 
0.3925059539229590 
0.4443099124772396 

10 0.0059745260011464 
0.0373611383696360 
0.0571010082270104 
0.1041446643790700 
0.1355161859798696 
0.2110062905707306 
0.2508943643743651 
0.3352721033248206 
0.3812593890562080 
0.4762946103276755 



exp(— ii/ir). exp{—iH2T) . exp(— 2ii73r). exp{—iH2T). exp(— ii/ir) (55) 
1 - 2ir (sBo + 5^tot.5tot) - 2r^ (sBq + Stot-^tot) ' • (56) 



Again, the effective Hamiltonian depends to second order only on the total spin 
component operators. The first and second order 5*3 DFS DD sequences can be expanded 
similarly and also decouple quantum baths. 

For third and higher order, both the ^3 and 5*3 sequences fail to decouple the 
quantum bath to the nominal order. To find higher order sequences for decoupling 
the quantum bath, we have resorted to brute force computational searches. We seek a 
sequence of Hamiltonians -ffcr(fe) and associated time intervals r^ such that the product 
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of unitary propagators 

N 

U = W exp(-iiJ^(fc)rfc) (57) 

fc=i 

has an effective Hamiltonian that couples only to the total spin operators Stot,xi Stot,y, 
and Stot,z, to the desired decoupling order. Here A^ is the total number of evolution 
intervals and cr{k) specifies under which Hamiltonian ( H3|) -(l48l) interval k evolves. 

Given some candidate sequence of Hamiltonians, we expand (157|) to the desired 
decoupling order in r. Because there are three qubits in the DFS, the system-bath 
interaction terms are at most weight three Pauli products on the system (constituent) 
qubits. The set of 64 Pauli products of weight three or less is partitioned into 
16 equivalence classes under the action of the group 5*3. For example, Xi is in 
the equivalence class consisting of [Xi,X2,X3], and X1Y2 is in the equivalence class 
consisting of [X1Y2, -^1^3, ^1-^2, ^1-^3,-^2^3, ^2-^3]- In the order-by-order expansion of 
fl57|) the coefficients of the Pauh products in the same equivalence class must be equal in 
order to achieve an effective Hamiltonian that couples to global interactions alone. Note 
that the coefficients consist of products of bath operators, with some numerical prefactor. 
Because the various bath operators do not commute, they further grade the terms that 
must be set equal. A candidate Hamiltonian sequence is a valid decoupling solution if 
Tfc's can be found such that members of each of the equivalence classes have the same 
coefficients, with the r^'s real and positive. For the examples above, if the members 
of the [^i]s3 equivalence class all couple to the bath operator C, the corresponding 
global interaction term is Stot,x-C- Similarly, if the members of the [XiF2]s'3 equivalence 
class all couple to the bath operator D, the corresponding global interaction terms are 

('-'tot,x'-'tot,i/ ~ I'^tot.z)--^- 

With this methodology, we have found the third order quantum bath decoupling 
sequence shown in table HI Rather than giving switching times, we have listed the 
evolution time intervals corresponding to a total evolution time of 1. The sequence 
consists of 26 Hamiltonian intervals, with a "doubly palindromic" structure. The first 
13 Hamiltonians are all even permutations while the second 13 are all odd permutations. 
The second 13 Hamiltonian types can be determined from the first 13 via the mapping 
Hi ^ H4, H2 ^ Hq, and H^ ^ Hr,. The interval timings for the odd permutations are 
identical to the even permutations. Within each set of 13 intervals, the interval timings 
are palindromic. Other length 26 third order solutions have been found, but the one 
displayed in table H] is optimal in the sense that the ratio between the maximum and 
minimum interval lengths is smallest. 

5. Numerical Simulations 

Numerical simulations of a DFS qubit coupled to classical and quantum baths have 
been performed to illustrate the performance of the DFS decouphng sequences. Figure 
[3] shows a log-log plot of the infidelity versus total evolution time for a DFS qubit coupled 
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Table 4. Third order DFS decoupling sequence for a quantum bath, consisting of 26 
Hamiltonian intervals. The first column lists the Hamiltonian type (H51) - (P5|) . The 
second column gives the interval length (rather than the switching times). The third 
column gives the pulse required at the end of the Hamiltonian interval. Only the first 
13 intervals are given; the interval timings and pulses for the second 13 intervals are 
identical. 



Hamiltonian type 


interval length 


pulse 


Hi 


0.02443154605193963 


P 


H2 


0.03273388118971666 


P 


H3 


0.05269740572865081 


p-i 


H2 


0.03073701555573789 


P-I 


Hi 


0.04633548169315730 


p-i 


Ha 


0.05049836419256131 


P 


Hi 


0.02513261117647280 


P 


H2 


0.05049836419256131 


p-i 


Hi 


0.04633548169315730 


P-I 


H3 


0.03073701555573789 


P-I 


H2 


0.05269740572865081 


P 


H3 


0.03273388118971666 


P 


Hi 


0.02443154605193963 


Pia 



to a classical dephasing only bath, with A3 DFS DD sequences of orders n = to n = 4 
applied to the DFS qubit. The infidelity is defined as 1 — F, where F is given by f lHT]) . 
The classical bath functions have size \Bj{t)\ ~ lOOMHz and \{dBj/dt)/Bj\ ~ lOOMHz. 
Each point in figure E] gives the infidelity after a single round of n**^ order decoupling 
for the corresponding total evolution time T on the abscissa; each point is averaged 
over 50 dephasing bath instances and 100 initial encoded DFS states. The plot markers 
correspond to A3 DFS DD orders n = through n = 4 as described in the figure 
caption. For n**^ order decouphng, the phase differences ( 134|) -( 136|) at short total times 
T have dependence ~ j^n+i^ ^]-^g corresponding infidelity has total time dependence 
~ 2"2(n+i)_ ^Yie dotted lines in figure [3] are fits to the function a — bT'^^"'^^^ for each order 
of decoupling. At short total times the infidelity scales as expected, with a straight- 
line dependence on the log-log plot. The slopes of the fits have the predicted values of 
2(n + 1), indicating that the DFS DD sequences decouple the dephasing baths to their 
designed orders. 

Figure m shows a log-log plot of infidelity versus total evolution time for a DFS qubit 
coupled to a quantum bath. We simulate a DFS qubit coupled to a bath composed of 
six spins; each constituent system qubit is coupled to two bath spins, and all six bath 
spins are coupled to each other. The Hamiltonian simulated is 

(2 4 6 \ 6 6 

j=l j=3 j=5 J j=l k=j+l 

which has the same structure as ( 142|) . In ( 158|) J gives the energy scale of the system-bath 
coupling, P gives the energy scale of the intra-bath coupling, Vij are random numbers 
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Figure 3. Simulations of a DFS qubit coupled to a classical dephasing-only bath. 
Infidelity is plotted against total evolution time T for ^3 DFS DD orders n = (free 
evolution) through n = 4. n = infidelities are given by • , n = 1 by ■, n = 2 by 
♦ ,n = 3byA,n = 4byT. Dotted lines are fits to the given leading order infidelity 
term. 
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between and 1, 5*^- is the vector of Pauli operators for the j*^ system spin, and Ij is the 
vector of Pauh operators for the j*^ bath spin. We take J = lOOMHz and /3 =10KHz 
for the simulations in figure HI The simulations are "numerically exact" in that the 
propagators for full system and bath evolution are determined by computing the matrix 
exponential of the Hamiltonian ( l58ll : swaps between system qubits are interspersed 
between the propagators as required by the DFS DD protocols. Fidelities are again 
computed according to (l3T]l . with the states, propagator U, and projector 11 suitably 
generalized to include both the system and bath, and the addition of a partial trace 
operation over the bath spins. Each data point in figure H] shows the infidelity computed 
for a given total evolution time, for decoupling orders n = through n = 3. Each data 
point gives the infidelity averaged over 52 random initial conditions and Hamiltonians. 
The first and second order decoupling sequences used are 5*3 DFS DD sequences; the 
third order sequence used is given in table HI As in the classical bath case, the infidelity 
is expected to scale as T'^(^~^'^\ The dotted lines giving fits to the function a — 5T^("+^^ 
confirm the expected scaling. 



6. Conclusion 



We have shown that exchange pulses alone are sufficient to decouple the 3-qubit DFS 
from classical and quantum baths. By averaging over permutations of the 3 constituent 
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Figure 4. Simulations of a DFS qubit coupled to a quantum bath. Infidelity is plotted 
against total evolution time T for DFS DD orders n = (free evolution) through n — ?>. 
n = \ and n — 2 use S3 DFS DD pulse sequences; n = 3 uses the pulse sequence in 
table m n — Q infidelities are given by • , n = 1 by ■, n = 2 by ♦, n = 3 by A. 
Dotted lines are fits to the given leading order infidelity term. 
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qubits, local baths can be made to appear global to high order. Because the 3-qubit 
DFS is immune to global decoherence, DFS DD protects the encoded information. 
Numerical simulations of the new DFS decoupling pulse sequences confirm DD sequence 
performance expected from analytical considerations. 

Decoupling of the 3-qubit DFS from classical and quantum baths may also be 
accomplished using NUDD (nested Uhrig dynamical decouphng) pulse sequences [TO] . 
The DFS DD pulse sequences, however, are far more efficient. For example, for third 
order decoupling from a quantum bath, n = 3, NUDD requires in + 1)2><3 _ 4ggg p^lse 
intervals, compared to 26 pulse intervals for the DFS decoupling described in section 
m Accounting for the structure of the DFS — protecting the encoded information only 
and using exchange pulses only — substantially reduces the number of pulses needed for 
decoupling and removes the need for single qubit Pauli gates. The fact that only a 
particular subspace or subsystem needs to be protected should be used in designing 
efficient decoupling pulse sequences for general qubit encodings. 

Decoupling by averaging over the symmetric group (or one of its subgroups) using 
exchange pulses can be generalized to other decoherence free subspaces and subsystems 
[I5] . The most efficient sequences as well as the types of errors protected against will 
depend on the encoding, the structure of the system-bath coupling, and the noise model. 
For the 2-qubit decoherence free subspace, for example, in which the encoding protects 
against collective decoherence in a single direction (say z), decoupling over 5*2 protects 
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against bath variations in the z direction only. Bath variations in the x and y directions 
will cause leakage out of the encoded subspace. Decoupling over S2 has exactly the 
same structure as UDD for a single qubit: a full swap on the exchange gate between 
the two qubits takes the role of the Pauli X pulse for the single qubit, and the interval 
timings are the UDD timings (E]). For other encodings protecting against weak collective 
decoherence [15], averaging over Sn subgroups will decouple bath variations in the 
encoding-protected direction only. For encodings protecting against strong collective 
decoherence (which includes the 3-qubit DFS) averaging over S'„ subgroups decouples 
the encoded information from all bath components. The effective global interaction 
created by the averaging affects only the gauge while leaving the encoded information 
unchanged. The quantum numbers describing encoded states correspond to total spin 
operators on increasing numbers of the constituent qubits, and all total spin operators 
commute with global interactions. Determination of switching times and pulse interval 
Hamiltonian types for decoupling other DFS encodings can be found by generalizing the 
methods described in this paper. 

Finally, we note that the correspondence between 5*2 decoupling and UDD of a 
2-level system, each coupled to classical dephasing bath(s) along a single direction, 
generalizes to a correspondence between decoupling a (i-qubit DFS and a rf-level system 
{SU{d)) from classical dephasing baths along a single direction [2l]. In these cases 
averaging over the d cyclic permutations of the system for the DFS or cyclic permutations 
of the d levels yields the same pulse interval Hamiltonian types and switching times. 
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Appendix 

For an initial encoded state \fpe) = (r, \/l — r^e'*^), the coefficients in the fidelity 
expression ( l32l) are 

Co = - (3 - 2r^ + 1r^ + 1r^{\ - r^) cos 20) , (A.l) 

6 



ci = -(1 - r^) [\ + 2r2 + 2rv/3(l - r^) cos^) , (A.2) 

C2 = -(1 - r^) {\ + 2r2 - 2r v/3(l - r'^) cos <\^ , (A.3) 

C3 = — (1 - 2r^ + lOr^ - 6r2(l - r^) cos 20) . (A.4) 
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